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Combining recent results on rational solutions of the Riccati-Schrodinger equations for shape 
invariant potentials to the finite difference Backlund algorithm and specific symmetries of the isotonic 
potential, we show that it is possible to generate the three infinite sets (LI, L2 and L3 families) of 
regular rational solvable extensions of this potential in a very direct and transparent way. 

PACS numbers: 



I. INTRODUCTION 



In quantum mechanics, the rarity of the potentials which are exactly solvable in closed-form (most of them belonging 
to the class of shape-invariant potentials^—) gives a undeniable importance to the reseach of new families of such 
potentials. A possible way to generate new solvable potentials is to start from the known ones and to construct regular 
rational extensions of them. If the procedure has a long history, in the last years important progress have been made 
in this direction^ - — . In a recent worki^ we proposed an approach allowing to generate such regular extensions 
' starting from every translationally shape-invariant potential (TSIP) of the second category (as denned in2£). For 
this, we use regularized excited states Riccati-Schrodinger (RS) functions as superpotentials in a generalized SUSY 
partnership. The regularization scheme corresponds to a "spatial Wick rotation" which eliminates the singularities 
from the real axis, a device already suggested by Shnol— in 1994 as a way to generate rational extensions of the 
harmonic potential. In the following years, this suggestion has been developped by Samsonov and Ovcharov^ and 
Tkachuls^ 3 -. Recently Fellows and Smith^i rediscovered this technique in the case of the harmonic oscillator, the second 
rational extension of which being the so-called CPRS potential^. Ini2, we have extended the procedure to cover the 
whole set of TSIP belonging to the second category. For the isotonic oscillator, we recovered the LI family of rational 
extensions discovered by Gomez-Ullate, Kamran and Milson^— , Quesne^r— and Odake, Sasaki et a l 13 ' 16 ' 17 . For the 
other second category potentials, the infinite set of regular quasi-rational extensions that we obtain coincides with 
^ i the Jl fiunilyftifc-iaiiil. 

In the present article, combining the finite difference Backlund algorithm with new regularization schemes which 
7—1 ' are based on specific symmetries of the isotonic potential, we show how the extension of the SUSY QM partnership 
t-H . to excited states allows to generate the three infinite sets LI, L2 and L3 of regular rationally solvable extensions of 
the isotonic potential (as well as the singular LO and L3 ones) in a direct and systematic way. This approach leads 
to a simple and transparent proof of the shape- invariance of the potentials of the LI and L2 series. 

The paper is organized as follows. We first recall how the generalization of the SUSY partnership based on excited 
states leads to a series of singular rational extensions of the initial potential. We then introduce basic elements 
concerning the finite difference Backlund algorithm viewed as a set of covariance transformations for the class of 
Riccati-Schrodinger equations and we interpret the generalized SUSY partnership in this perspective. In the third 
and fourth sections, we recapitulate some results concerning the isotonic oscillator, its connection with confluent 
hypergeometric equation and the Kienast-Lawton-Hahn's Theorem which describes the distribution of the zeros of 
the Laguerre functions on the real axis. The fifth section is devoted to present the set of parameters transforms which 
are discrete symmetries of the isotonic potential. Using them as regularization transformations, we show then that 
the finite difference Backlund algorithm based on the corresponding regularized RS functions generates directly the 
three series LI, L2 and L3 of regular rationally solvable extensions of the isotonic potential. In the last section, we 
prove the shape-invariance of the potentials of the LI and L2 series. 



II. GENERALIZED SUSY PARTNERSHIP BASED ON EXCITED STATES: LO SERIES OF RATIONAL 

EXTENSIONS 

Consider a family of closed form exactly solvable hamiltonians H(a) = —cP/dx 2 + V(x; a), a £ W n , i e J C 1, the 
associated bound states spectrum of which being given by (E n (a), w n (x; a)), where w n (x;a) — — ip' n (x; a)/ip n (x; a) 



2 



is the Riccati-Schrodinger (RS) function associated to the n th bound state eigenfunction ip n (x;a). The Riccati- 
Schrodinger (RS) equation^ for the level E n (a) is then 

- w' n {x; a) + wl(x; a) = V(x; a) - E n (a), (1) 

where we suppose Eo(a) = 0. The RS function presents n real singularities associated to the n simple nodes of the 
eigenstates ip n ( x 'i a )- As it is well know n 26 ' 27 , H(a) admits infinitely many different factorizations of the form 

H(a)-E n (a)=A+(w n )A(w n ), (2) 

where 

A(w n ) = d/dx + w n (x;a), (3) 

with, in particular 

A{w n )^ n (x;a) =0. (4) 
This allows to associate to H(a) or V(x; a) an infinite family of partners given by 

# (rl) (a) - E n (a) = A (w n ) A+ (w n ) = -d 2 /dx 2 + V {n) (x; a), (5) 

with 

V {n) (x;a) = V(x;a) +2w' n (x;a). (6) 

For n > 1, these potentials are all singular at the nodes of tp n (x; a) and are defined on open intervals only. 
On these domains, H^ n \a) is (quasi)isospectral to H(a). Indeed, writing 

^\x;a) = A(w n )ip k {x;a), (7) 

it is easy to verify that we have for any k 

H^\a)^\x;a)=E k {a)^\x;a), (8) 
that is, ip^ l \x;a) is an eigenstate of H^ n \a) associated to the eigenvalue Eu(a). We write symbolically 

V (n) {x:a) ee V(x;a), (9) 

iso 

where = means "isospectral to" . Defining 

iso 

w n , k (x;a) = - ^ X;a) (10) 



Eq.(|SJ| gives, for k > n 



v' nM {x;a) +w n , fe (x;a) 2 = V (n) (x; a) - E k (a). (11) 



This scheme generalizes the SUSY QM partnership, by using the excited state RS functions w n as superpotentials. 
However, only for the ground state n — 0, the factorization and then the partner potential (a) = V(x; a)+2w' (x; a) 
are non singular and we recover the usual SUSY QM partnership^. 
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III. FINITE DIFFERENCE BACKLUND ALGORITHM 

We can consider the preceding partnership in a different way which gives a prominent role to the covariance 
transform of the RS equations class. 

A. Invariance group of the Riccati equations 

As established by Carihena et aL 28 i 29 , the finite-difference Backlund algorithm is a consequence of the invariance of 

the set of Riccati equations under a subset of the group Q of smooth SL (2, R)- valued curves Map(R, SX(2,R)). For 
any element A £ Q characterized by the matrix: 

A{x) = ( jU S(x) ) ' det A{x) = a{x)S{x) " /?(x)7(a;) = X ' (12) 
the action of A on Map(R, R) is given by: 

n A~ n a{x)w{x) + 
w(x) — > w(x) = —— — =7 — r-. (13) 

If A acts on a solution of the Riccati equation: 

w'(x) = ao(x) + ai(x)w(x) + a2{x)w 2 (x) , (14) 
we obtain a solution of a new Riccati equation: 

w'(x) =ao(x)+ai(x)w(x)+a2(x)w 2 (x), (15) 
the coefficients of which being given by 



u 2 {x) 
ui(x) 
u (x) 



~a (x) = M{A)~t{x) + $(x), lt(x)=\ ui(as) I, (16) 



where: 

/ 5 2 {x) -l{x)S{x) 7 2 W \ / W(j,6;x) 

M(A) = -2/3(x)d(x) a(x)S(x) + j3(x)-y(x) -2a(x)j(x) , W(x) = W(6, a; x) + W{/3, 7; x) \ (17) 

\ (i 2 {x) -a(x)P(x) a 2 (x) J \ W(a,/3;x) 

(W(f,g;x) — f(x)g'(x) — f'(x)g(x) is the wronskian of f(x) and g(x) in x). As noted inM, Eq. (fl6l) defines an affine 
action of Q on the set of general Riccati equations. 

B. Particular case of the RS equations and finite difference Backlund algorithm 

The most general elements of Q preserving the subset of RS equations has been determined in 2 ^. Among them we 
find in particular the elements of the form: 

AM -^(«? »><>, (18) 

where cj)(x) satisfies an RS equation with the same potential as in Eq.([T]) but with a shifted energy: 
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~<j>'{x)+cf{x) = V(x)-(E-X). (19) 
With this choice w(x) satisfies the RS equation: 

- w'(x) + w 2 (x) = %{x) - A, (20) 

where %{x) = V(x) + 2<f>'(x). 

Consequently, starting from a given RS function of the discrete spectrum w n (x;a), for every value of k such that 
E k > E n , we can build an element A (w n ) € Q of the form: 

A( \ 1 ( w n (x;a) E k (a) - E n {a) -w n 2 {x;a)\ 

H> y/E k (a)-E n (a) V - 1 w ^ a ) J 

which transforms w k as: 

, ,A{w) ( „) E k (a)-E n (a) 

w k (x;a) -> u>l '(a; a) = —w n {x; a) H - -, (22) 

w n (x\a) - w k (x;a) 

where w k is a solution of the RS equation: 

- w k n) \x-a) + (w k n) (x-a)y = V^(x;a) - E k {a), (23) 
with the same energy Ei(a) as in Eq([T]) but with a modified potential 

V {n \x;a) = V(x;a) + 2w' n (x;a). (24) 

This is the content of the finite-difference Backlund algorithm^— . It transposes at the level of the RS equations 
the covariance of the set of Schrodinger equations under Darboux transformations^—. In the following we call A (w n ) 
a Darboux-Backlund Transformation (DBT). 

To V(x; a), A (w n ) associates the (quasi)isospectral partner V^(x; a). Among the A (w n ), only A (wq) leads to the 
regular, usual SUSY QM partner V^(x; a). The correspondence between the eigenvalues of V(x; a) and V( n \x;a) is 
direct. We also have from Eq(j7|) and Eq((3|) 

ip { k n \x; a) ~ (w n (x;a) - w k (x; a)) ip k (x; a), (25) 



w n ,k (x;a) = w k (x;a). (26) 



that is (see EqljlO 



Then, the finite difference Backlund algorithm generates exactly the RS functions corresponding to the spectrum 
of the generalized SUSY partner V^ n > of V . 

Note that for shape invariant potentials (SIP)**, A (wq) is in fact an invariance transformation of the RS equations 
associated to the considered family of potentials (indexed by the multiparameter a), since in this case 

V {a) {x;a) = V(x;ax) + R(a) (27) 

and 

w { k °\x,a) = w k _i(x,ai), (28) 
where a\ = /(a) and R(a) axe two given functions of the multiparameter a. 
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As we noted before, starting from the RS function w n of a regular excited bound state which has n nodes on the real 
domain of definition I of V(x;a), we generate via A(w n ) a generalized SUSY partner which presents n singularities 
on this domain. Nevertheless, the finite difference Backlund algorithm can be applied by replacing w n by any other 
solution of the same RS equation Eq([T]), even if this solution does not correspond to a physical state. Knowing 
w n {x, a), the general solution of Eq([T]) is given by 

2 / x x w n (s;a)ds 

W n (x;a,W ) = w n (x;a) - & - Wn{t:a)dt , (29) 

Wo + J Xo dse ™ v ' 

where Wo is an arbitrary real parameter. We could then use the DBT A (W n ) to build a generalized SUSY partner 
potential V^ n \x; a, Wo) = V(x; a) + 2W' n [x\ °j W>) and look for values of Wq for which W n and are not singular. 
For some potentials it is nevertheless possible, by using specific symmetries, to build directly the researched regular 
RS functions. Such symmetries exist in particular for the isotonic oscillator. 

IV. THE ISOTONIC OSCILLATOR 

As shown in 20 -, the primary translationally shape invariant potentials (TSIP), for which a\ = a + a, can be classified 
into two categories in which the potential can be brought into a harmonic or isotonic form respectively, using a change 
of variables which satisfies a constant coefficient Riccati equation. 

The first clement of the second category is the isotonic oscillator potential itself (ie the radial effective potential for 
a three dimensional isotropic harmonic oscillator with zero ground-state energy) defined on the positive real half line 

V(x;u,a) = ^-x 2 + a ^ ? ^ +V (u,a), x > 0, (30) 
4 x z 

with a = I + 1 > 1 and Vq(uj, a) = — uj (a + j)- The shape invariance property of V(x; ui, a) is expressible as 

U (0) (x;u;,a) = V (x; lo , ax) + 2uj (31) 

and its spectrum is given by 



E n (uj) — 2nu), ip n (x; u>, a) ~ exp y — J w n (x; u>, a) dxj , (32) 

where the excited state Riccati-Schrodinger function (RS function) w n (x;w,a) can be written as a terminating con- 
tinued fraction as 



with 



w n (x; u), a) — wq(x; uj, a) + R„ (x; u>, a), 



i \ w a 
wo(x;uj,a) — — x 

2 x 



(33) 
(34) 



and 



R n (x;uj,a) = - 



E n (uj) 



w (x;uj,a) +w n -i(x;uj,a 1 ) 

—2nuj 2 (n — j + l)uj 



(35) 



2u) 



ujx — (2a + 1) jx- 



ujx - (2 (a + j) - 1) fl- 



ux — (2 (a + n) — 1) fx 



As it is well known, the isotonic oscillator eigenstates can be also expressed in terms of Generalized Laguerre 
Polynomials (GLP) as 
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<p n (x; w, a) - /e- ul!/4 4"- 1/2) (wa; 2 /2) . (36) 

This implies that we have 

R n {x^a) = - (log (4°-^ (- 2 /2)))' = ^ ^/^r^l ( 37 ) 

L„ (fjJX / Z) 

which is singular at the nodes of ip n (x;u>, a), that is, at the zeros of L)? (£). Concerning these last, we have a 
classical result of Kienast,Lawton and Hahr*2Zr— : 

Kienast-Lawton-Hahn's Theorem 

Suppose that a £ — N. Then Ln (z) admits 

1) ft positive zeros if a > — 1 

2) n + [a] + 1 positive zeros if — n < a < —1 ([\oe\] means the integer part of a) 

3) No positive zero if a < — n 

The number of negative zeros is always or 1 . 

1) 0ifa>-1 

2) if -2k - 1 < a < -2k and 1 if -2k < a < -2k + 1, with -n < a < -1 

3) if n is even and 1 if n is odd, with a < —n 

This theorem, confirms in particular that, for positive values of a, the RS function w n (x;uj,a) corresponding to a 
physical bound state and then the associated generalized SUSY partner V^ n '(x;uj,a) of V(x;uj, a) present always n 
singularities on the positive half axis. This family of singular rational extensions of V(x;u),a) will be called the LO 
series. 

V. CONFLUENT HYPERGEOMETRIC EQUATION AND ISOTONIC OSCILLATOR 

The confluent hypergeometric equation 

zy (z; a, A) + (a + 1 — z)y (z; a, A) + Xy(z; a, A) = (38) 

on the positive half real line, can always been reduced to a Schrodinger equation for an isotonic oscillator. Indeed, if 
we put z = ujx 2 /2 and <fi (x; a, A) = y(z; a, A) in Eq. (|38j) . we obtain the following equation for (f> (x; a, A): 

<f>"{x; a, A) + ( 2 " + 1 - wx)<f>'(x; a, A) - 2uX(f>{x; a, A) = (39) 



Then 



satisfies 



il>(x; a, A) ~ <\> (x; a, A) exp(i f dxfi^- - ujx)) = x a+112 e"^ x ' !/4 (x; a, A) (40) 

2 J X 



- i/}"(x; a, A) + ( + ( Q + 1 /2)(a — l J2± _ u t a + u\ M X] a \\ = (x; a, A) . (41) 
\ 4 x z J 

If we define a = a + 1/2, ip(x; a — 1/2, A) = ipx( x > a ) an d E\(lu) — 2Xlu, we obtain 

H(u, a)<il>x(x; a) = £ A (w)^ A (x; a), (42) 

where 

i/. A (x;a) ~ x a e-^ 2 / 4 ?;(a;a; 2 /2; a - 1/2, A), (43) 
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H(uj,a) being the usual isotonic hamiltonian (see Eg. (130])) 

H(u,a) = -^ + V(x;u,a). (44) 

Eq. (|4"2")) is the Schrodinger equation for the isotonic oscillator, where for physical bound states we must have A = n. 
In this case, the confluent hypergeometric equation 

zy"(z]a- 1/2, n) + (a + 1/2 - z)y'(z; a - 1/2, n) + ny(z; a - 1/2, n) = 0, (45) 
admits the regular solution 

y(z- 1 a-l/2,n) = Li a - 1 / 2 \z) (46) 

and we have 

ip n (x;a) ~ x a e-^ x2 l i L^- 1 ' 2 \iox 2 /2). (47) 
This is exactly the physical state for the isotonic oscillator at the energy E n — 2nuj. 

In fact, as shown by Erdely i 4 ' 8 ' 38 ! 39 , Eq . f|38]) admits quasi rational solutions built from GLP in four sectors of the 
values of the parameters a and A 



\ = n: y (z;a,n) — L ( n\z) 



A = — n — a — 1 : yi(z; a, a + 1 + n) = e z Ln (—z) 

A = n — a 
X = -n-l: y 3 (z;a,-n- 1) 



y 2 (z;a,n - a) = z a Ln a \z) 



(48) 



*e*Lk~ a \-z) 



They correspond to the following four eigenfunctions 



i> n (x;a) ~ x a e- ux2 [ / A L [ n~ 1,2) (ujx 2 /2), E n (u) = 2nuj 
*p- n - a -i/2(x; a) ~ x a ef** 3 /*L£- 1/2 \-ua?/2), E_ n _ a _ 1/2 (u;) = -2 (n + a + 1/2) u 
Tp n -a+i/2(x; a) ~ x^e'"* 2 /*Ln la ~ 1/2) (wx 2 /2) ) E n _ a+1/2 {u) = 2(n-a + 1/2) u 
V>-„-i(x;a) - x 1 - a e" x2 / 4 L,T (a_1/2) (-wa; 2 /2), £L„_i(w) = -2(n+l)w 



(49) 



The 3 last cases don't correspond to physical states and physical energies. 



VI. DISCRETE SYMMETRIES OF THE ISOTONIC RS EQUATION 



Since the isotonic oscillator is shape invariant, the A (wq) DBT is an invariance transformation for the RS equations 
associated to the family of isotonic oscillators indexed by the couple of parameters (w, a) (see Eq. (|30| ). But this family 
of RS equations is covariant under other specific transformations which act on the parameters of the isotonic potentials 
and preserve their functional class. As we will see, the connections between the quasi rational sectors of the confluent 
hypergeometric equation admit a very simple interpretation in terms of covariance transformations of the isotonic 
potential. 



A. Inversion of ui the parameter 



The first covariance transformation for V(x;u>,a) acts on the u) parameter as 



W F4(_ W ) J V(x;u,a)^V(x;u,a)+u(2a+l) 

1 w n (x;u),a) -4 v n (x;uj,a) = w n (x;-u),a), 
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v n (x; uj, a) satisfying (E n (— a;) = —E n (uj) = E_ n (uj)) 

- v' n (x; uj, a) + v 2 n (x; uj, a) = V(x; uj, a) - £L( n+0+1 / 2 ) (uj) . (51) 
From Eq.®, Eq.flM]) and Eq.©, writing 

v n (x; uj, a) = v Q (x; uj, a) + Q n {x; uj, a), (52) 

we deduce 

. . uj a 

v (x;ui,a) = -—x (53) 

2 x 

and 

Qn{x; U ,a) = — (54) 

v (x;w,a) + v„_i(x;w,ai) 

2nuj 2 (n — j + 1) uj 2uj 



+ (2a + 1) /x+ "' wi + (2 (a + j) - 1) /x+ " + (2 (a + n) - 1) jx 

log(^-V 2 ) (- W x 2 /2)))'. 



Clearly, for a > 1 (Z > 0), v n {x;u>,a) does not present any singularity on the positive real half line. This result is 
coherent with the above mentioned Kienast-Lawton-Hahn's theorem since the argument of the GLP Ln in the 
expression of Q n is now a strictly negative value. 

Note that we recover exactly the same results if we use the " spatial Wick rotation" 19 i 21 ~ ^i 

w n (x; uj, a) — > v n (x; uj, a) — iw n (ix; uj, a). (55) 

This means that the T u transformation send the singularities of w n , which are initially all on the real axis, on the 
imaginary axis. This explains why the new RS function v n does not present any singularity on the real line. Finally, 
comparing Eq. (|4TJ]) to Eq. (|5"2")) . Eq. (|53"|) and Eg . (j5"4")) . we see that transforms an eigenfunction of the first sector into 
an eigenfunction of the second sector and then coincides with the Kummer's transformation^. 

B. Inversion of a the parameter 

The second covariance transformation acts on the a parameter as 

a *4 1 - a, { ^ Ow) ^ V (x; w, a) + u(2a - 1) (5g) 



w n (x;ui,a) -4 u n (x;uj,a) = w n (x;uj, 1 - a), 



u n (x;ui,a) satisfying 



- u' n (x; uj, a) + u 2 n (x; uj, a) = V(x; uj, a) - E n+1/2 _ a (uj) . (57) 
From Eq.Q, Eq.flM]) and Eq.([35]) we deduce 



u n (x; uj, a) = u n (x; uj, a) + P n (x; uj, a), (58) 

where 

u (x;uj,1) = %x+ - — - (59) 
2 x 

and 
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Pn{x; U ,a) = — , (60) 

v (x; ui,a) + v n -i(x;uj,a-i) 

—2nuj 2 (n — j + 1) u) 2lj 

f ... f A- J ; ' , f* ... I* 



ujx + (2a — 3) fx— ujx + (2 (a — j) — 1) jx— ujx + (2 (a — n) — 1) /x 



If in this case the argument of the GLP in the right hand member is strictly positive, the associated a = — (a — 1/2) 
parameter being strictly negative. In accordance with the Kienast-Lawton-Hahn's theorem, by taking a sufficently 
large, we can decrease the number of real zeros and in particular we can eliminate all the positive zeros. Thus, if 
a > n + 1/2, L n t " a 1 (ujx 2 /2^ is strictly positive for any value of x. This means that P n (x;oj,a) and u n (x;u),a) arc 
not singular on ]0, +oo[ when a = n + m + 1/2, with m > 0. 

Note that 

P 1 (x;uj,a) = — —— = Q 1 (x;uj,a-2). (61) 

ujx + (Za — 6) /X 

Finally, comparing Eq. (l49l) to Eq. (l58l) , Eq. (l59l) and Eq. (1601 , we see that T a transforms an eigenfunction of the first 
sector into an eigenfunction of the third sector. 

C. Inversion of both parameters u> and a 

Finally, we can also act simultaneously on both parameter as 

( W , a) F ^ l-a)\ V ^ ^ V ^ ^ a) + 2UJ (62) 

{ w n (x;u),a) a -s>" r n (x;u),a) = w n (x; -u, 1 - a), 

r n (x;u>, a) satisfying 

- r' n (x;u,a) + r 2 n (x;uj,a) = V(x;-u, 1 - a) - E n (-u) = V(x;uj,a) - E_ (n+1) (uj). (63) 
From Eq.([3"3"|). Eq.® and Eq.([55|) we have 

r n (x;u),a) = r (x;uj,a) + T n (x; u>, a), (64) 

where 

. . w a — 1 . 

ro(x; lj, a) = x H = — w (x; u), a — 1) (65) 

2 x 

and 



T n (x;uj,a) = — En ^ (66) 

r (x;uj,a) + r n _i(x; lj, a_i) 

2nu> 2 (n — j + 1) u> 2lo 



wi — (2a — 3) /x+ uix — (2 (a — j) — 1) /x+ ujx — (2 (a — n) — 1) /x 

-(lo g (47 (a - 1/2) (-^ 2 /2)))'. 



In this case, the argument of the GLP in the right hand member and the associated a = — (a — 1/2) parameter are 
both strictly negative. In accordance with the Kienast-Lawton-Hahn's theorem, by taking a sufficently large, we can 
have any zero on the negative half line if n is even (n = 21) and one if n is odd. Thus, if n — 21 and a > 21 + 1/2, 
L 2 ^ a (— lox 1 12) is strictly positive for any value of x. This means that Tziix; co, a) and r 2 ;(x; co, a) are not singular 
on ]0, +oo[ when a = 21 + m + 1/2, with m > 0. 
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Note that 



2(jj 

Ti(x;u,a) = — r- = Rx(x\w, a - 2). (67) 

ujx — (2a — 6) fx 



Finally, comparing Eq. (|49| to Eq. (|64[) , Eq. (|65|) and Eq. (|66| , we see that r a o T u transforms an eigenfunction of the 
first sector into an eigenfunction of the fourth sector and corresponds also to a Kummer's transformation^. 

VII. REGULAR RATIONAL EXTENSIONS OF THE ISOTONIC OSCILLATOR 

Since the transformations considered above are covariance transformations for the family of isotonic potentials which 
regularize the RS functions, we can use these regularized RS functions into the finite difference Backlund algorithm 
and generate regular isospectral partners for the isotonic potential. 

A. Rational extension of the LI series 

Wk and v n are associated to the same potential but with different eigenvalues (cf Eq (|5T|l ) 

-v' n (x;ui,a)+vl(x;uj, a) = V(x;uj,a) - E_ (n+a+1/2) (uj) , gg . 



-w' k (x;uj, a) + w\(x; uj,a) — V(x:uj,a) — E% (uj) . 
which means that we can use v n to build a DBT A (v n ) and apply it to w k as 

, {n), \ , x E k (u>) ~ £L( n+a +i/ 2 )(w) 

w k (x;uj,a) -> wl (x; uj, a) = —v n (x; uj, a) H -. r -. -, 

v n (x;uj,a) - Wk(x;uj,a) 

where w k n \x;uj,a) satisfies 



(69) 



- w k n) ' (x; u,a)+ (w k n) (x; uj, a)) * = V {n) (x; uJ,a)-E k (uj) , (70) 

with 

V (n) (x;uj,a) = V(x;uj,a) + 2v' n (x;LO,a). (71) 
For every n > 0, (x; uj, a) is regular on the positive half line and isospectral to V(x; uj, a) 

V {n) (x;uj,a) = V(x;uj,a). (72) 

iso 

Clearly, vJl\x; uj, a) = —v n (x;uj,a) is also a solution of Eq (|70|) associated to the eigenvalue E_/ n+a+ i/ 2 )( ljJ ) < 

Eq (uj) = 0. Nevertheless, its asymptotic behaviour is similar to the one of w_ (x; uj, a) — —lux/2 — a/x and conse- 
quently 

ip^ (x; uj, a) ~ exp ( — / w^l\x; uj, a)dx ) (73) 



cannot satisfy the boundary condition associated to the physically allowed eigenstates. 

All the physical eigenfunctions of H^ n > (uj,a) = —d 2 /dx 2 + V^ n ' (x; uj, a) are then of the form 

ip k n) (x;uj,a) = r—^. = . A(v n )tp k (x;uj,a), fc>0 (74) 

y/E k (uj) - E_( n+a+1 / 2 )(uj) 
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and is strictly isospectral to H. 
Since (cf Eq®) 

V(x; oj, a) + 2v' (x; oj, a) = V(x; oj, at), (75) 
Eq(|7T|) and Eg ([72]) can still be written as 

V {n) {x;oj,a) = V(x;u>,ai) + 2Q' n (x; oj , a) = V{x;oj,a). (76) 

iso 

For instance, we have for n = 1 

T^fiv \ t r / \ 4w 8w(2a+l) 

wa^ + 2a+l (ojx 2 + 2a + l) 

and we recover the first rationally-extended radial oscillator obtained by Quesnei^. For n = 2, we have immediately 
from Eq. (J52J, Eq.flSJ) and Eq. dMJ) 



cj a-1 4wa; (wi 2 + (2a + 1)) 

2 X x (uj X 2 + ( 2a +l)) 2 -2(2a+ 1) : 



w 2 (x;cj,a- 1) = -xH h - — — — — ——2 — — — — , (78) 



which corresponds to the superpotential associated to the second rationally-extended radial oscillator of the LI series 
obtained by Quesne^. 
More generally, we have 

Q n (x;oj,a) = (log (l^ 1 ^ (-ojx 2 /2)))' . (79) 
In Odake-Sasaki 's approac h 13 ' 16 ' 17 , this corresponds to a prepotential of the form 

W n (x;oj,a) = -^x 2 + alogx + log (l^ 1 / 2 ) (-ojx 2 /2)^j (80) 

and we recover (up to a shift ino->a + ii-2) the result obtained i n 13 ' 16 ' 17 ancU^ for the potentials associated to the 
LI exceptional orthogonal polynomials. 

B. Rational extension of the L2 series 

As in the preceding case (cf Eq ([571) h we can use u n to build a DBT A (u n ) 

( {n) ( \ / E k {oj) - E n+1/2 _ a {oj) 

w k (x;oj,a) -> w k [x; oj, a) — —u n {x; oj, a) H - 7 -, (81) 

u n [x;oj, a) - w k (x;oj,a) 

where w k (x;oj,a) satisfies 

- w ( k n) '{x- oj, a) + (w k n) (x; oj, a)) * = U^(x; oj, a) - E k (oj), (82) 

with 

U {n \x;oj,a) = V(x;oj,a) + 2u' n (x;oj,a). (83) 
If a > n + 1/2, (x; oj, a) is regular on the positive half line and isospectral to V(x; oj, a) 
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U (n \x;uj,a) = V(x;w,a). (84) 

iso 

In this case, as for the LI series, we see immediately that w_ (x; uj, a) = —u n (x; uj, a) is another solution of Eq(|82|) 



associated to the eigenvalue E n+ x/2- a (oj)(uj) < Eq (uj) = 0. But here again, the asymptotic behaviour of w_ (x; uj, a) 
is similar to the one of w^(x;uj,a) — —ujx/2 — (a — l)/x and consequently 



ifj_ '(x;w,a) ~ exp I — J w_ (x;ui,a)dxj (85) 

cannot satisfy the boundary condition associated to the physically acceptable eigenstates. 
All the physical eigenfunctions of H^ n >(uj, a) = — d 2 /dx 2 + U^(x;uj,a) are then of the form 

ip^\x;uj,a) = r — ! -. = ===A(u n )tfj k (x;uj,a), k>0 (86) 
^E k (uj) - E n+ i/ 2 _ a (u) 

and in the L2 series, H^ n ' is also strictly isospectral to H. 
Since (cf Eq. ((59])) 



V(x; uj, a) + 2u' Q (x; uj, a) = V(x; uj, a_i), (87) 

using Eq.([83|) and Eq.([84|), we obtain 

U {n) (x;ui,a) = V(x;u),a-i) + 2P^(x;u,a) = V{x;ui,a). (88) 

iso 

Note that, since P\(x;uj,a) — Q\(x;ui,a — 2), the first rational extension of this family has the same functional 
form than the first rational extension of the preceding family. 
For instance, we have for n = 2 

Aujx (ujx 2 + (2a - 5)) , , 

P 2 (x;uj,a) = i \ '1 , (89) 

+ (2o-5)) + 2(2a-5) 

which corresponds to QuesneAS second rational extension of the L2 series. 
We have also, by redefining a — > n + a 

V(x; u, a„_i) = V(x;uj,a n ) + 2i^(i;w,a„), (90) 

iso 

where 

p n (x-u,a n ) = — — ; J — ? ... r — — — — - (91) 

ujx + (2n + 2a — 3) fx— ujx + 2 ((n + a — j) — 1) fx— ujx + [2a — 1) fx 



log(^(«+»- 1 /2) (ujx 2 12)))' 



is regular on the positive half line for a > 0. In Sasaki and a l 14 i 17 formulation, we recover the associated prepotential 
via 

W n (x; u,a) = - J u n (x; uj,a + n)dx = -^x 2 - - + " - 1 + log (£-(«+»-V2) (ujx 2 /2)^ (92) 
and the family of regular rational extensions obtained is exactly the L2 one. 
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C. Rational extension of the L3 series 

Finally, Wk and r n being also associated to the same potential but with different eigenvalues (cf Eq (|63[) ). here again 
we can use r n to build a DBT A (r n ) and apply it to ioj. 

, ^fcO („) E k (uj) - E_ (n+1) (uj) 

w k (x;uj 7 a) -t wl J (x;uj,a) = —r n {x;w,a) H r, (93) 

r n (x;w,a) - Wk(x;u,a) 

where w k (x; u>, a) satisfies 

- w { k n) '(x; u, a) + (w[ n) (x; u>, a)) * = (x; w, a) - E k (u), (94) 

with 

W {n \x;uj,a) = V(x;io,a) +2r' n (x;uj,a). (95) 
If n — 21 and a > 21 + 1/2, W( 2 *) (a:; w, a) is regular on the positive half line and isospectral to V(x; u), a) 

W {2l \x;oj,a) = V(x;u,a). (96) 

iso 

As for the eigenfunctions of H (21 \lo,o) = -d 2 /dx 2 + (x; u, a) generated from those of by the DBT Eq ([93]) . 
they are given by 

^ l \x;w,a)= 1 ==A(r 2l )il> k {x;u>,a), fc>0 (97) 

y/E k {w) - E_ {2l+1) (w) 

and constitute physically allowed eigenstates. But for the L3 series the isospectrality is no more strict as for the 
preceding series. Indeed, Eg ([Ml) is evidently satisfied by the regular RS function 

w_(x;u),a) = — r2i(x',cj,a), (98) 
the asymptotic behaviour of which being identical to the one of (x; w, a) = — ro{x;u),a). Then 

and contrarily to the preceding cases 

ip^ l \x;u},a) = exp (J dxw^L l \x;cj,a)\ (100) 

is a physical state associated to the eigenvalue E_^ n+1 ^(uj) < 0, that is, the fundamental state of the hamiltonian 
H^(ui,a). Consequently, H^ 2V> and H are only quasi-isospectral in this series, admitting a supplementary 

energy level lower than those of H. 
Since (cf Eq.(|g5])) 

V(x; co, a) + 2r' (x; uj, a) = V(x; to, a_i) — 2a;, (101) 
using Eq. ([95|) and Eq. (|96|) . we obtain 

W (n) (x;uj,a) = V{x]Lo,a_ l )-2Lo + 2P' ri {x- 1 Lo,a) = V(x;u,a). (102) 
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Since T\(x; lo, a) = Ri(x; lo, a — 2), the first rational extension of this family has the same functional form than the 
first rational extension of the LO family. 
For n = 2, we have 

-4lox (lox 2 - (2a - 3)) , , 

T 2 (x;co,a) = i ^ ^ , 103 

(lux 2 - (2a-3)) 2 + 2(2a-3) 

which is regular if a > 2 (I > 1) and corresponds to Quesne 1 ^ second rational extension of the L3 series. 
If we redefine a — > 21 + 1/2 + a, 

T 21 {x;lo,21 + cl+1/2) = - log fe/° +2i) (-wx 2 /2))' (104) 

and W-^ 2 ') (x; w, a + 21 + 1/2) are regular on the positive half line for a > 0. 

VIII. SHAPE INVARIANCE PROPERTIES OF THE EXTENSIONS OF THE ISOTONIC OSCILLATOR 

As observed initially by Quesne&iS' on the n = 1 and n — 2 examples, the rational extended potentials of the LI and 
L2 series inherit of the shape invariance properties of the isotonic potential. Several general proofs of this result have 
been recently proposed 1 -^^, in particular by Gomez-Ullate et al£. In the present approach, these shape invariance 
properties can be derived in a very direct and transparent manner. 

A. Shape invariance of the extended potentials of the LI series 

The superpartner of a potential of the LI series VW{x;u,a) = V(x;lo,o) + 2v' n (x; lo, a) is defined as 

V (n) {x;Lo,a) = V {n) {x;uo,a) + 2w { q 1) ' {x;lo,o), n > 0, (105) 

Wq 1 \x\ lo, a) (see Eq. (p9"|) ) being the RS function associated to the ground level of (E Q (lo) = 0). 
We then have 

V*0(*; W ,a) = yW(x;,,a)-2<(x;,,a)-2( E_ {n++1/2) (co) V 

\v n (x;to,a) - w (x; lo, a) J 



V(x;lo,o) - 2 



E -(n+a+l/2)(u) 



v n (x;uj,a) - wo(x;LO,a) t 
Using Eq(|53]). the shape invariance property of V(x; lo, a) in Eq. ([31]) can also be formulated as 

V(x; lo, a) + 2v' Q (x; lo, a) — V(x; lo, a%). (107) 
Inserting Eq (fT07)) in Eq (fT06|) . we obtain 

V^\x;Lo,a) = V(x-Lo, ai )-2( g -(»+°+V2)H + ^ a) \ (10g) 

\v n (x;uj,a) - w (x;Lo,a) ) 

= V^(x;Lo,a 1 )~2(Al l )' , 

where 



E -(n+a+l/2){u) . 

1 vo{x; lo, a) + v n (x; lo, a\). (109) 



n v n (x;uj,a) - w (x;Lo,a) 
As an example, consider the special case n — 1. Using Eq (|54p . we can write 
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A{ = -2u (a + 3/2) * ujx- + — , (110) 

vq(x;uj ,a)-\-Vo[x',U} ,eti) 

ujx+^ 2a +1 2ujx 
= (2a + 6) — 5 — — - — uix ~ — — = —ujx. 

y 1 ujx 1 + (2a + 3) x ujx 2 + (2a + 3) 

We obtain 

V (1) (x;uj,a) = V {1 \x;uj, ai ) +2uj, (111) 

which implies that V^(x;uj,a) has the same shape invariance properties as V(x;uj,a). 
More generally, using Eq(|54]) and defining z — —ujx 2 /2 and a = a + 1/2, we obtain 

A n = E -(a+ n +\/2)(u) / . 1 n _ + {vo(x;u),a) + v (x;uj,a 1 )) + Q n (x\uj,a x ) (112) 

= 2a + 2 ifr-Vjz) (z) + Lj a) (z) 2a 

x L i n a \(z) + L^- 1) (z) l£\z) 

But the generalized Laguerre polynomials satisfy the identity 

L^(z) + L^(z) = L^(z), (113) 



which gives 



A i _ _ (a + n) ifr-V (z) + zL { n a+1) (z) - qljg (z) 

zL\ (z) 

The other fundamental recurrence 

(n + a) L { :\ (z) - zL^+V {z) _ ( „ _ z) L M {z) = 0> (115) 
combined with Eq (|113[) gives then directly 

^ = -wi, (116) 

that is, 

V {n) (x;uj,a) = V ( - n) (x;uj,a 1 ) + 2uj. (117) 
Consequently V^ n '(x;u},a) inherits of the shape invariance properties of V(x; uj,a) for every value of n. 

B. Shape invariance of the extended potentials of the L2 series 

The superpartner of a potential U^ n '(x; uj, a) — V(x; uj, a) + 2u' n (x; uj, a) of the L2 series is defined as 

U {n) (x;uj,a) = U (n) (x;uj,a) + 2w { n) '(x;uj,a), n>0, (118) 
Wq U \x; uj, a) (see Eq. (|5Tj) ) being the RS function associated to the ground level of . Then 
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UW(x;u>,a) = V(x;u>,a)-2( g " + . 1/2 " fl( f r) . (119) 

Using as before, the shape invariance properties of V(x; ui, a), this gives 

tT (n) (!6iU>o) = V(x-u: 1 a 1 )~2( ^n+i/2-aM .^u^ (120 ) 

= D-W(x; W ,oi)-2(Aj)', 

where 

A n =- B n+l/2-o(w) 7 \~~ 7 7 + Vq(x; OJ, a) + U n (x] OJ, Oi). (121) 

Using Ea(|54p and defining z = ujx 2 /2 and a = 1/2 — a, this becomes 

(2n-2a+l)w / L^ (z) (n + a)L ( n a) (z) 



= r~ — 4 — r + Pn(x: oj, at) — ujx I . " 1 , v ' — I ^— — ' — - — , v I , (122) 

But the generalized Laguerre polynomials satisfy the identity 

zl£+» (z) = (n + a)Li a \ (z) nL™ (z) , (123) 

which combined to Eq (|113p gives 

Al = -ujx. (124) 
Then U {n \x;uj,a) has the same shape invariance properties as V(x; uj, a) for every value of n, that is 

U (n) (x; ui, a) = U [n) (x; u), o x ) + 2w. (125) 

C. SUSY partners of the L3 series extended potentials 

In this case, the superpartner of the extended potential V^(x;uj,a) = V(x;u,a) + 2r' n (x; oj , a) is defined as 

W {n \x;uj,a) = W {n) {x;uj,a) + 2(-r' n (x;u,a)) = V{x;uj,a), n>l (126) 

since —r n (x; ui, a) is the RS function associated to the ground level of W^ n \ 

The SUSY partner of W^ n \x;uj,a) is nothing but the initial potential V(x;uj,a) itself and the DBT A(v n ) is the 
reciprocal of a SUSY partnership. 

IX. CONCLUSION AND PERSPECTIVES 

In this article, a new method to generate the regular rational extensions of the isotonic oscillator associated to the 
LI and L2 families of exceptional Laguerre polynomials is presented. It is based on first order Darboux-Backlund 
Transformations which are built from excited states RS functions regularized by using specific symmetries of the 
isotonic potential. Starting from this primary shape invariant potential and using the combination of these symmetries 
and DBT (as covariance transformations), we generate four towers of secondary potentials, the four series LO, LI, L2 
and L3. Among them, the potentials belonging to the LI and L2 series are regular as well as half of the potentials 
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of the L3 series, the other ones being singular on the positive half line. The secondary potentials of the LI and L2 
series inherit of the same translational shape invariance properties as the primary isotonic potential. 

These new potentials being obtained, it is still possible to use the Krein-Adler theore m 40 ! 41 and its subsequent 
extension obtained by Samsonov 42 , to generate other secondary potentials by applications of some particular n th 
order DBT. 

A similar study can be conducted for the other second category potentials (Darboux-Poschl- Teller or Scarf hyperbolic 
and trigonometric) but also for the first category potentials. These laslj^ - include the well known case of the one- 
dimensional harmonic oscillato r 4 ! 19 ! 21- — but also the Morse potential (the regular algebraic deformations of which 
having already be obtained by Gomez- Ullate et al 4 -), the effective radial Kepler- Coulomb potential and the Rosen- 
Morse potentials. This work is in progress and will be the object of a forthcoming paper. 
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